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ABSTRACT 

A simple model for the hadronic contribution to the photon vacuum polarization 
function YlhadiQ"^), for spacelike momenta, is presented. For small momenta, the two 
loop contribution from the pseudoscalar meson octet is computed from the chiral La- 
grangian. The light quark contribution (which at low momentum gives the 0{q^) 
counterterm in the chiral Lagrangian) is calculated within a relativistic constituent 
quark model incorporating the momentum dependence of the quark mass. The per- 
turbative gluons of QCD are included in a standard fashion. The total result is close 
to an estimate of HhadiQ"^) that is obtained directly from e~^e~ — *• hadrons data. We 
further use our results for lihadiq^) to calculate the C(e^) hadronic contribution to 
lepton magnetic moments and to calculate aQED{,Mz). A simpler model of constituent 
quarks with momentum independent masses gives less favourable results. 
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Introduction 



The photon vacuum polarization function, n(g^), where 

m^^(g) = ie'^Ii{q^)[q^g^y - q^q^] (1) 

plays an important role in any high precision QED observable, e.g. the lepton anoma- 
lous magnetic moments ai {I = e,fi,T), atomic energy levels, e'^e~ — >• e~^e~ scattering. 
For many practical applications, one requires n(g^) for spacelike momenta only (g^ < 0) 
so we will specialize to this region, where 11 (g^) is a smooth and real function. We define 
our subtraction point such that n(0) = 0. 

The largest source of uncertainty lies in the hadronic contribution, Uhadil'^)- The 
most reliable estimates of IlhadiQ'^) that are available use experimental data for 

a{e~^e~ — 7* — hadrons) 
^ a{e+e- ^ 7* ^ /^V") ^ ' 

in conjunction with dispersion relations. A simple analytic fit to this estimate is given 
in lU, and we will compare our result to this fit. 

From the theoretical viewpoint, UhadiQ"^) is in principle calculable in terms of the 
(current) quark mass parameters and a QCD scale (e.g. A;^). Our approach is a 
less ambitious one, namely to extract UhadiQ"^) (and other related physical observables) 
from a simple QCD-inspired model. We will show that a simple picture can go a long 
way towards providing some theoretical understanding of Uhadiq'^)- 

The behavior of HhadiQ"^) at both low and high momenta is constrained by QCD. 
The high momentum description is provided by perturbative QCD. The low momentum 
dependence receives an important contribution from pseudoscalar meson loops, and 
this can be calculated using standard techniques of chiral Lagrangians. The relevant 
quantity occurs at 0{q^) in the low energy derivative expansion, and thus an associated 
two- loop calculation must be performed. Knowledge of an 0{q^) counterterm is also 
necessary, and at present any estimate of this quantity is model dependent. 

In this paper we shall be using a relativistic constituent quark model of the three 
light quarks to smoothly interpolate between the low and high energy QCD contribu- 
tions to Uhadil'^)- This is a gauged nonlocal constituent (GNC) quark model which 
incorporates the momentum dependence of the quark mass as a natural regulator. The 
pion decay constant and all other quantities appearing in the chiral Lagrangian to 
0{q^) have been expressed in terms of this mass function. In particular the standard 
quantities 0] Li, L2, Lg, and Liq are well described in terms of one parameter (de- 
noted by A below) appearing in the mass function|0]. Reasonable values for the other 
Lj's and the current quark masses are obtained as welljQ. The model has also been 
successfully applied to certain other quantities, the pion electromagnetic form factor 
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and the vector-minus-axial two point function, beyond There has also been 

a recent discussion [|| of the advantages of our nonlocal regularization in the context of 
anomalous processes. These results suggest that some of the more essential aspects of 
nonperturbative QCD are accounted for by the momentum dependence of the quark 
mass. 

There are two GNC quark loop diagrams which contribute to UhadiQ"^)- At high 
momenta the effects of the constituent quark mass become negligible and the GNC di- 
agrams smoothly approach the naive one-quark-loop contribution to UhadiQ"^)- We add 
perturbative QCD corrections to the model description of the three light quarks. More 
precisely we add the hard gluonic corrections beginning at 0{as) as extracted from 
a standard perturbative QCD calculation of Uhadil'^)- The contributions to UhadiQ"^) 
from the c and b quarks will be represented completely by perturbative QCD. 

We will consider the implications that our calculated UhadiQ"^) has for the hadronic 
contributions to the muon and tau magnetic moments. Of the two, the muon magnetic 
moment is less affected by the various uncertainties. We will also consider the QED 
running coupling at large momentum, aQEniMz). 

Before discussing the various components of our calculation in more detail, we 
present the main results in figs. |l] and ^ The effects of the meson loops and the 
hard gluonic corrections are shown separately. They are added to the GNC quark 
contribution to produce our total result for UhadiQ"^)- We find that the total result is 
quite close to the experimentally-based estimate]^]. We note as well that for small — g^, 
the GNC quark-loop contribution is significantly larger than a more naive constituent 
quark model with momentum independent masses, for example with the typical values 
rriu = rrid = 330 MeV and = 550 MeV. This comparison is shown in fig. ^. 

Although we divide our calculation into various parts, we stress that we are con- 
sistently using one model to describe Uhadil"^) for all spacelike momenta. This is the 
GNC model with perturbative gluonic effects added. The GNC model includes the 
pseudoscalar mesons, and as we have said, it nicely reproduces the standard chiral 
Lagrangian of low energy QCD. We will therefore describe our calculation of Uhadil"^) 
at low — using the language of chiral Lagrangians. In the model all meson dynam- 
ics, including the meson kinetic terms, are generated through quark loops. At higher 
energies the compositeness of these mesons will become evident via form factors, and 
their further contribution to UhadiQ"^) will be damped out. We will enforce this feature 
of the model by cutting off the meson-loop contribution at a scale of order rrip. 
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Pseudoscalar mesons 



For —q^ much less than mp, Iihad{<f) can be evaluated from the SU (3) l x SU{?>)r chiral 
Lagrangian with explicit symmetry breaking terms. The one-meson-loop contribution 
corresponds to a calculation at 0{q^) in the low energy expansion and the counterterm 
at this order is eliminated by the definition Ilhadi^) = 0. The leading behavior of 
^had{(f) requires a calculation at 0{q^), and this corresponds to a two-loop calculation. 
It is this piece which contains the leading logarithms ln{fi/mTr) and ln{fi/mK)- We will 
add this two-loop piece to the finite one-loop contribution, and the result should provide 
a good approximation for Ilhad{(l'^) at low — g^. All diagrams with more than two meson 
loops are of 0{q^) or higher, and are therefore neglected. 

We will follow the definition of dimensional regularization counterterms [LKfi)) 
given in |0. We do not enforce the 0{q^) equations of motion, so we must reinstate 
two terms 0] that are usually removed by the equations of motion. In the end our result 
can be expressed in terms of the L[(/x)'s of as well as one new parameter, C"'(/i), 
which is the 0{q^) counterterm. 



lihad{q^) 

n(4) ^.2^ 



had 




dz In 



mt 



dz In \ 1 — zil — z) 



m 



K , 



had 



3(47r/o) 



+ ri(/i) + T2 + Ts{q') + T,iq' 



(3) 

(4) 
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/o is the pseudoscalar decay constant in the chiral limit, and explicit expressions for 
the Tj's are given in the appendix. We require that the dependence of C^{fJ^) cancels 
that of Ti(yu). 

To obtain the numerical value of Ilhad{(l'^) for small — g^, we need only the values 



of the coefficients C"'(/i) and L'[{^). We will choose /i = mp and take the experimental 



imp)'s from 



]. We will give below the value of 



values of the L\' 
quark model. 

But to isolate a purely meson- loop contribution to Ilhad{ 



[nip) from the GNC 



we will remove the 



nir 



from d^). This term will effectively be included in the low momentum behav- 



ior of the quark-loop graphs of the GNC quark model. The remaining terms in (^) 
represent the two-loop meson contribution. We will allow this plus the one-loop result 
to contribute to the growth of IlhadiQ'^) up to —q^ = m^. On these scales there may be 
substantial error in a two loop calculation; we consider the implications of this below. 
The one and two loop meson contributions are displayed in fig. and compared to 
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the quark-loop contribution. The latter dominates even though the contributions at 
low — are inversely proportional to the mass squares of the quark and pion. (The 
relative enhancement of the quark loops is due mainly to the colour and spin degrees of 
freedom.) But the meson- loop contribution is still significant. We find an uncertainty 
of 20% in the meson-loop contribution at due to the uncertainty in the L['s. 

The fact that the one-loop meson contribution is smaller than the two-loop contri- 
bution is consistent with the dominance of the p meson. The p is, of course, integrated 
out of the chiral Lagrangian which means that its effects are incorporated into the 
counterterms L[, C", etc. Since these counterterms do not appear in (H), the one- loop 
meson contribution is suppressed. 



GNC quarks 

The GNC Lagrangian contains the octet of pseudoscalar mesons vr" and a quark 
triplet ip with a dynamical quark mass and a current quark mass matrix M.. 
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Note that X(x, y) is a path-ordered exponential. = — ^^75 and i?^ = + ^^75 
are left and right handed external gauge fields, respectively. For = the model has 
local SU{?i)L X SU{?i)R symmetry, hke QCD in the presence of external gauge fields. 
Numerically, we use the current quark masses m„ = = 8 MeV and = 180 MeV. 

The dynamical quark mass S(— g^) in (0) is the Fourier transform of the S(x — y) 
appearing in (0). The parameter A specifies the value of mo through its relation to /q. 
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For the most part we will use the values /o 
mo = 317 MeV. 



84 MeV and A = 2 which correspond to 
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The GNC quark contribution comes from two one-quark-loop contributions to the 
vacuum polarization — one with the two external photons attached at one vertex, 
and the other with the two photons attached at two distinct vertices. The resulting 
contribution to HhadiQ^) is shown in fig. In this contribution we have included the 
small effects of the naive one-loop graphs of the c and b quarks. (The commonly-used 
ranges 1.3 < rric < 1.7 GeV and 4.7 < < 5.3 GeV produce an uncertainty less than 
2%.) The gluonic corrections to light and heavy quark-loop graphs will be treated 
below. We may also consider the sensitivity of the GNC result to the quantities /o 
and A within an allowed range 84 MeV < /o < 88 MeV and 2 < A < 3.|1 The GNC 
contribution is reduced by less than 15% and we find that 0.19 < C""(mp) < 0.20 . 

In fig. ^, we compare our m, d, s GNC result to the analogous contribution from a 
more naive model with momentum independent masses of typical values m„ = = 
330 MeV and = 550 MeV. Note that the shapes of the two curves are quite different, 
and this is true for any values of the momentum independent quark masses. The GNC 
curve is more consistent with the fit of Burkhardt et. al.|jl| 



Pert ur bat ive gluons 



The final contribution that must be considered is due to radiative gluons. We have al- 
ready included one-quark-loop diagrams with the correct non-zero masses (momentum 
dependent masses for the light quarks), so we wish to extract the gluonic corrections 
given by terms containing a^. We use the following dispersion relation 

= Zf- /■°°^,^(^)-^o(s 

pert 



127r2 



s s 



where q < 0. R{s) has been calculated using the MS scheme to 0{{as/TT) ) in 
Nf quark fiavours with charges Qi. 
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Ro{s) is the value of R{s) when is set to zero, and its appearance in 
the naive one-quark- loop result. 

as{s) is obtained by solving the QCD P function^. 
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(3o = 11 -^iV/ (17) 

A = 102 -^Nf (18) 

2857 5033 ,^ 325 ,. 

02 = Nf^ Nf (19) 

2 18 ^ 54 ^ ^ ' 

We have cut off the integral in ([T^) at s = . The result is an estimate of hard gluonic 
corrections missing in the GNC model, and it is shown separately in figs. |^ and 0. The 
gluonic contribution contains an uncertainty of roughly 23%(16%) at — = m^(M|) 
corresponding to the range .11 < as(M|) < .13 and an uncertainty of 26%(6%) at 
-g2 = m2(M|) for 1.3 < < 1.7 GeV. 

By our choice of the cutoff in (|l^) we have added to the GNC quark model only those 
perturbative corrections which can be reliably calculated. Our intent is to see how the 
model does at describing the contributions that cannot be calculated perturbatively. 



Lepton anomalous magnetic moments 

To lowest order in uqed the hadronic contribution to a lepton's magnetic moment has 



the following algebraic form 110 
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The superscript "(4)" reminds us that this is all at 0{e^). Notice that < for the 
entire range of integration, so we only need the value of HhadiQ^) for spacelike g^. 

We perform the double integration by fitting our numerical result for Uhadil"^) to 
a piecewise-analytic function of < 0. Since -C rric, our estimated corrections 
to IlhadiQ'^) have essentially no effect on a^^^^^^^. For the case of a^^j^^^ the uncertainties 
in the meson-loop contribution and the corrections are both larger, with the former 
dominating. We take a 40% error for the meson-loop contribution in this case, double 
the naive error due to the uncertainty in the L['s as noted above, to include possible 
corrections from higher order effects in the chiral Lagrangian. The results are 



fi,had 

(4) 
^T,had 



(6.3 ±0.5) X 10"^ 
(3.2 ±0.1) X 10"*^ 



(23) 
(24) 
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Any other choice of {fo,A) in the range (84 MeV < /o < 88 MeV, 2 < A < 3) would 

reduce our resuh for a^'^Ld(^T^Ld) by less than 10%(6%). 

These results should be compared to the calculation from integrating an experimen- 
tal determination of R{q^). {R{q^) is defined in (Q). The precise form of the integral is 
given by Kinoshita et. al.|lTT|) 



f ajfi . = (7.05 ± 0.08) X 10-8 [11] 
experiment =^ I fi^'^ ^ ^ , LJ 25) 

' <L = (3.6 ± 0.3) X 10-6 111 ^ ' 



2\ 



On the other hand the parametrization of the experimentally-determined HhadiQ 
due to Burkhardt et. al.[|ll may be used directly in eqs. (E0|-B3). This gives 



afi„rf = 6.63 X 10-8 
experimental fit ^ <^ fi^'^ , 26) 

^ - 3.45 X 10-6 



''T,had 



The authors of n seem to claim that the uncertainty should be less than 5%, making 
this result for a^^^ noticeably smaller than the preceding result of ([25|). 

Finally, we consider quarks with momentum independent masses of = = 330 
MeV and = 550 MeV, and add the pseudoscalar mesons, heavy quarks and 
corrections exactly as discussed above. 



rrid = 330MeV, m, = 550MeV I 41'^'^ j^'^ ^ ^'^J ' (27) 



ii'L = (2.2 ± 0.1) X 10-6 



Clearly the GNC model is a significant improvement. In order to reproduce the GNC 
results for both lepton magnetic moments in (|23|) we would require different sets of 



masses, for example niu^^s = 243 MeV and rriu^^s = 201 MeV for a^'^jj^^ and a^^j^^^ re- 
spectively. This illustrates the fact mentioned previously that momentum independent 
quark masses cannot reproduce the shape of our Iihad{,(f) for the entire range —q^ < 0. 



QED running coupling 

It is straightforward to evaluate the hadronic contribution to the running of aQED{~q^) 
from our results. |TB| Using /o = 84 MeV and A = 2 we find 
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47r[n 
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Tlhadm = 3.68 ± 0.07 (28) 



had 
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The error is dominated by the error in the perturbative QCD contribution. This may 
be compared to the result in 

A = 3.94 ±0.12 (29) 

When combined with the well-known non-hadronic effects, our result ( pHD implies 

«qed(M|) = 129.05 ±0.07 (30) 

We note that this is also consistent with an independent determination using recent 
LEP data. (In that analysis additional corrections are included to define an effective 
coupling, aeff{Mz)). Any other choice of {fo,A) in the range (84 MeV < /o < 88 MeV, 
2 < y4 < 3) would reduce the value of A by less than 2%. 



Comments 

In this paper we have used a nonlocal constituent quark model for the description of 
^hadil^) at spacelike momenta. The success of our quark level description could be 
considered to be a manifestation of "duality" , and it clearly relies on being far removed 
from the resonance structure of QCD appearing for timelike momenta. The various 
uncertainties we have quoted occur within the model itself, and they are not intended 
as an a priori estimate of how well the model should resemble QCD. It is only after a 
comparison with experimental data that our model is able to shed light on some of the 
essential features of QCD dynamics. 

After completion of this work, we received a preprint [^] containing an indepen- 
dent theoretical estimate of the hadronic vacuum polarization and the muon magnetic 
moment. 



Appendix 

The complete result for YlhadiQ"^) for small — can be derived from the chiral La- 
grangian, and is given in (|]). Here we provide the explicit form of the Tj parameters, 
using the notation of at the renormalization scale fi. (The fi and dependence of 
all parameters is implicit in this appendix to simplify the notation.) The parameters 
rn-TT, and are the lowest order mass values. The relation of these parameters to 
the physical masses can be taken from eq. (10.7) of 0. The Tj's are well-behaved func- 
tions of spacelike such that g^Tj vanishes when g^ = 0. The apparent dependence 
in T2, T3, and T4 cancels out. 

Ti(/i) = 8l^ + 8lK-:r7^{ll + ll + ljK) (31) 
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Figure Captions 



Figure 1: All contributions to our total result for UhadiQ"^) for < SOOMeV. The fit 
from HH is also shown. 



Figure 2: All contributions to our total result for Ilhad{(l'^) for |g| < 3GeV. The fit from 
|l|] is also shown. 



Figure 3: Comparison of the contribution to HhadiQ"^) from the u, d and s quarks in 
two distinct approaches: the GNC model (as used in the present analysis) and a simple 
model with momentum independent masses m„ = = 330 MeV and = 550 MeV. 
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This figure "fig2-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9304264v2 
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